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PRESSURE MODEL OF A FOUR-WAY SPOOL VALVE FOR SIMULATING 
ElfCTROHYDRAULIC CONTROL SYSTEMS 
by Vernon D. Gebben 
Lewis Research Center 

SUMMARY 

Two programs for modeling the spool- valve portion of electrohydraulic servovalves 
are described. One computer model represents the valve as a flow device where spool 
position and load reaction pressure are input signals. It uses the underlap-spool-valve 
equation to define leakage through spool clearances. This flow model is a conventional 
representation for applications where fluid compressibility is the immediate system 
element loading the valve. The other computer model simulates a pressure source by 
using the opposite causality; that is, pressure is the dependent variable and flow is an 
independent variable determined by load reactions. It uses a new equation that provides 
the pressure-flow characteristics for spool leakage. This pressure model was designed 
for use in simulations that include hydraulic line dynamics. These models are de- 
scribed and defined by their equations, computer flow diagrams, and FORTRAN expres- 
sions. They have demonstrated characteristics that are realistic for time-domain 
simulation of electrohydraulic control systems. 


INTRODUCTION 

Electrohydraulic servovalves made their appearance in the late 1940*s to satisfy 
aerospace needs for compact, high-power, fast-response control systems. Analyses 
of these valves have remained involved and tedious. Difficulties are evidenced by the 
fact that articles on analysis techniques for estimating d5mamic performance of hy- 
draulic control systems are still being published in technical journals. 

Problems in valve simulation result from inherent nonlinearities. For example, 
the quadratic pressure-flow relation of orifices presents one problem source. Analysis 
problems also arise from nonlinear characteristics caused by spool overlap, leakages 
through spool clearances, and limited spool displacement. 


Analyses of hydraulic control systems have been structured around standard tech- 
niques of linear approximations and describing functions. Linear approximations ob- 
tained from the small-signal theory (theory of small perturbations) have provided solu- 
tions that proved particularity effective in evaluating system stability in the frequency 
domain (refs. 1 and 2). Stability studies in the frequency domain for large signals also 
have been accomplished by using describing functions that simplify the nonlinear prob- 
lem by excluding harmonic components from signals (refs. 3 to 6). Although these 
techniques have demonstrated successful applications, they become inaccurate when the 
system contains multiple nonUnearities such as a combination of coulomb friction and 
backlash. 

Analog and digital computer simulations are more direct and flexible. Computer 
models can represent most nonlinear electrohydraulic control systems with sufficient 
accuracy for time-domain and frequency-domain design analyses and provide the con- 
venience of permitting direct changes to parameters. 

For one example, a computer simulation of a hydraulic control system was chosen 
for evaluating d 5 mamic effects of long lines connecting the valve to the actuator. Ser- 
vovalves usually are simulated with a "flow model" that relates output flow as a function 
of both spool position and output pressure, where output pressure is determined by load 
reactions. The flow model for this application produced a program that used differen- 
tiators to simulate line dynamics. Differentiators, however, should be avoided when- 
ever possible, especially in analog circuits, where they often saturate computer ele- 
ments with amplified circuit noise. The "pressure model" described in this report 
provides an advantage by allowing line dsmamics to be implemented with integrators 
instead of differentiators. It is a model whose output pressure depends on the combina- 
tion of spool position and output flow. 

A pressure model that includes leakage through spool clearances was not available 
in the literature. Flow models seem to be the only kind reported. Furthermore, the 
nonlinear equations of the flow model could not be transformed into a pressure model. 

It was therefore necessary to invent a pressure model that closely approximates the 
pressure-flow characteristics produced by the conventional flow model. The new com- 
puter model is described in this report. 

Mathematical models for the spool portion of electrohydraulic servovalves are the 
subject of this report. The first section describes the pressure-flow characteristics 
of four-way spool valves and gives the flow-model equation. The second demonstrates 
the need for a new model. The new pressure model is then presented and compared 
with the conventional flow model. In the last section, the models are used to derive 
differences between the output characteristics of a servovalve connected to an inertial 
load and those of one connected to a capacitive load. Both models are defined in the 
appendix by their equations, computer flow diagrams, and FORTRAN expressions for 
digital programs such as the "Continuous System Modeling Program" (ref. 7). 
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SYMBOLS 


A amplitude of spool motion, dimensionless 

B blocked- line pressure gain (8p/9x) 

C normalized hydraulic capacitance or compliance 

Cj hydraulic capacitance or compliance 

D computer program variable 

E electrical potential 

F computer program variable 

G(s) transfer function 

H computer program variable 

I electrical cxirrent 

j imaginary unit, 

K computer program variable 

L dimensionless displacement parameter 

Mj inertia or mass 

N no-load flow gain (9q/9x) 

P normalized differential output pressure, dimensionless 
p pressure 

Q normalized output flow, dimensionless 

q flow 

R(X) modulated resistance, dimensionless 

Rj resistance 

s Laplace operator 

t time 

U normalized spool center region, dimensionless 

V normalized spool overlap, dimensionless 

W computer program variable 

X normalized spool position, dimensionless 

X spool position 
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y computer program variable 

Z computer program variable 

io angular frequency 


ELECTROHYDRAULIC SERVOVALVES 

A typical servovalve is shown in figvire 1. An input electrical- power signal oper- 
ates a torque motor. The torque motor operates a flapper valve which drives the spool 
to position X. Spool position controls the output pressure-flow relation. The main 
background information for this report is given in reference 8, which describes four- 
way valves, the pressure-flow characteristics of four-way valves, and the frequency- 
response characteristics of electrohydraulic servovalves. 

For analysis, the operating characteristics of servovalves can be separated into 
dynamic and static components as shown in figure 2 . The first block represents the 
combined dynamics of the torque motor, flapper valve, and spool. Generally, the dy- 
namic portion is assumed to be a linear third- or fourth- order system. The linear 
frequency-response characteristics that define G(s) can be obtained from the manufac- 
turer of the servovalve. In certain applications, it may be necessary to include flapper- 
valve nonlinearities which can significantly influence the d 3 ?namlc performance capa- 
bilities of fast-response systems. This limitation has been investigated and formulated 
into general design criteria (ref. 9). 

Function R(X) is the subject of this report. It is a load-dependent nonlinear func- 
tion comprised of orifice equations that relate the output pressure-flow characteristics 
to input spool position. Except for the overlap case, the present analysis assumes a 
symmetrical spool with port areas proportional to X. Also, as is commonly done to 
obtain a simplified analytical representation of R(X), it will be assumed that the flow 
in one output line always equals the flow in the other output line. Equal flows require 
identical d 3 mamic loads on the two output ports, which is a conservative estimate for 
stability analysis (ref. 10). The output flow can then be represented by a single vari- 
able, Q. 

Symbols Q, P, and X denote normalized flow, normalized pressure, and nor- 
malized spool position. Variable Q is equal to the average value of the flows through 
the two output ports divided by the maximum output flow at zero- differential output 
pressure. Variable P is equal to the differential output pressure divided by the max- 
imum differential output pressure that the valve produces for a blocked output Une. 
Variable X is normalized spool position defined as a fraction of maximum displace- 
ment. 
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Ideal critical- center (critical- lapped) four- way spool valves have characteristics 
that are displayed in fig^e 3 and defined by 

Q = X Vl - P if 0 < X < 1 (1) 

Q = X Vl + P if -1 X s 0 (2) 

The characteristics at null (X = 0) are very important to basic performance. Fig- 
ure 4 shows that the no-load (P = 0) flow sensitivity is linear and has a slope value of 1. 
Figure 5 shows the blocked- line (Q = 0) pressure sensitivity. Pressure P is either 
+1 or -1, which means that it is impossible to maintain zero pressure across the load. 
The ideal critical- center valve would thus be impractical for use in position-control 
systems containing coulomb friction. 

Practical spools have radial clearances and can have minute underlap allowed by 
machining tolerances. The result is bypass leakage for spool positions where X is 
less than position U. A typical value for U would be 0. 03; that is, 3 percent of max- 
imum spool displacement. The center region defined by U is small, but important to 
accvirate simulations. A small amount of leakage along the spool is beneficial, for it 
reduces blocked- line pressure sensitivity and provides intermediate output pressures, 
as shown in figure 6. Spools can also be manufactured with overlap. In this case, the 
pressure-flow characteristics at null will be a combination of a dead zone (threshold) 
caused by overlap and bypass leakage through radial clearances. 

Figure 7 has the pressure-flow characteristics divided into three zones. The ideal 
critical- center spool-valve equations (eqs. (1) and (2)) fit the principal regions very 
well. The center region generally has characteristics similar to those of underlap 
(open-center) four-way spool valves. The constitutive relation for the center region 
would then be described by the following relation, derived from reference 8, section 
5-4; 

Q = (0.5)(X+ U) Vl - P + (0.5)(X - U) Vl + P if -U s X s U (3) 

Equations (1), (2), and (3) are limited to a range of values of P between 1 and -1. 
Inertial loads, however, can produce pressure transients that put P outside of this 
range (ref. 11). High pressvme peaks occur when, for eJKimple, a valve switches from 
fully open (X = 1) to closed (X = 0) in an effort to stop a fast-moving mass. High differ- 
ential pressure then develops across the actuator as the result of high resistance of the 
closed valve and high flow maintained by load momentum. Since P is not physically 
limited between 1 and -1, a new set of relations between X, P, and Q are needed for 
simulating system transient response. 

Pressure limits in eqtiations (1) to (3) result from incorporating the following ori- 


5 



fice equation: 

q = D ypj - pg (4) 

where 
q flow 

D conductance coefficient 
Pj upstream pressixre 
Pg downstream pressure 

This equation provides no information concerning flow direction, except that upstream 
pressure is higher than downstream pressure. Conventional network equations of the 
four-way spool valve have fixed directions of flows (ref. 8, p. 79). 

A different set of equations describing the constitutive relations of four-way valves 
results when flow polarity is added to the orifice network equations (ref. 12); namely, 

Q = [sign (1 - P)](X) yil - P| if U < X s 1 (5) 

Q = [sign (1 - P)](0. 5)(X + U) y |l - P| + [sign (1 + P)](0. 5)(X - U) 

if -U < X ^ U (6) 

Q = [sign (1 + P)](X) y|l + P| if -1 < X s -U (7) 

Equations (5), (6), and (7) can be combined and expressed by a single constitutive 
relation: 

Q = [sign (1 - P)](0. 5)(X + L) yfl - P| + [sign (1 + P)](0. 5)(X - L) y|l + P| (8) 

where 

L = U if |x| ^ U 
L = |X| if |X| s U 

Figure 8 shows the pressure-flow characteristics obtained from equation (8) for 
pressures in the extended range where P is greater than 1. The scales of figure 8 are 
noticeably different from those of figure 6. To make the curves of equation (8) visible 
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in the center region, the graph was stretched in the Q-axis direction by multiplying 
variable Q by constant 1/U. Lines of X, designating spool position, were changed to 
Unes of X/U. 

Although these equations extend the pressure range, they are considered valid only 
if cavitation does not occiur. This restraint comes from the definition of Q which 
assumes that the two output ports have equal flows. The output of models described in 
this report should therefore be limited to conditions where P does not exceed 1. Equa- 
tion (8) still provides the advantage of allowing simulations to operate in the extended 
pressure range (P > 1) without program interrupts. The effects of excessively high 
values of P can then be qualified (ref. 12). 

Equation (8) is a flow model, designated ’'model Ql. " The appendix gives a com- 
puter program for equation (8) that provides the same press\ire-flow characteristics 
used in other published computer programs (refs. 11 to 15). Also included in the pro- 
gram is a dead zone that simulates the characteristics of overlap spools. A program 
that solves for Q, however, is not always best for system simulation, as will be 
demonstrated in the next section. 


MODEL OPTIONS 

Load elements of typical actuation systems shown in figure 9 include line inertia, 
fluid compliance, actuator leakage, actuator inertia, and friction. Not all of these 
elements, however, may be necessary for an accurate system analysis. For example, 
line d 5 mamics generally can be neglected if the valve is mounted on the actuator. Fluid 
compliance can also be neglected if the resonant frequency due to compliance is much 
higher than the system bandwidth. It will be demonstrated that the type of valve model 
needed for a simulation depends on the immediate element loading the valve. 

Dynamics of capacitive and inertial load elements can be expressed mathematically 
in one of two forms: either integral causality - 



q dt 


p dt 


(9) 

( 10 ) 


or derivative causality - 
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( 11 ) 


q = Cj 

p =Mj 


dt 

dt 


( 12 ) 


Interdependence of valve model and load equations can be demonstrated with linear 
load elements and a linear valve model. Let the valve be defined by the following equa- 
tion: 


X - P - 5L = 0 (13) 

B N 

where B is the blocked -line pressure gain, and N is the no-load flow gain. The 
causal form of the valve flow model is written 

q = N^x-£j (14) 

For the valve pressure model, 

p-B^x-^j (15) 

Two different programs result when these equations are applied to the load shown in 
figure 9. Figure 10 describes the simulation with the valve flow model, and figure 11 
describes the simulation with the valve pressure model. The simulation of figure 10 
requires a differentiator (M^s) to compute the pressure that accelerates the fluid in the 
line. Algebraic manipulation of figure 10 will not eliminate the differentiator. Differ- 
entiation disappears only when line inertia is zero. Then fluid compliance is the im- 
mediate element loading the valve. 

The simulation of figure 11 has the advantage of using an integrator (1/Mjs) to com- 
pute line flow. This causal form is preferred over that of figure 10, because computers 
are more accurate at integration than differentiation. Also, digital computer programs 
are difficult to stabilize when a set of simultaneous equations (algebraic loop) contains 
differentiation. Therefore, the valve pressure model should be used in applications 
where the nearest load element to the valve is inertia. 

Likewise, the valve flow model is preferred when line inertia is not included in the 
simulation. Actuator pressure, determined by fluid accumulated in the cylinder, is then 
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related to the time integration of the difference between valve flow and flow due to ac- 
tuator motion. 

Valve models with flow for their dependent variable seem to be the only kind re- 
ported. The author is unaware of any previous descriptions of pressure models (i.e. , 
valve models that use output pressure as the dependent variable). Hence, a pressure 
model representing four-way valves was developed for analyzing systems involving line 
dynamics and other systems where flow models are unusable. 


VALVE MODEL PI 

The pressure model defines P as a function of X and Q. The direct method for 
deriving this causal form would be to express the flow model in terms of P. Unfor- 
tunately, this approach is obstructed by square-root terms in equation (8). Equation (8) 
cannot be solved explicitly in terms of P. Therefore, a new constitutive relation ap- 
proximating equation (8) is needed for modeling four- way spool valves. 

To provide realistic valve characteristics, the new relation should have the follow- 
ing features: 

(a) Identical equations for the principal operating region ( | X | ^ U) 

(b) Pressure gain (3P/9X) that does not change abruptly during transition from prin- 

cipal operating region to center region 

(c) Realistic slopes (9P/9Q, 9P/9X, 9Q/9X) around the origin (X = 0, P = 0, Q = 0) 

(d) Equation that is tangent to equation (8) at the origin 

These conditions were satisfied by an equation that was invented, rather than having 
been derived from a set of theoretical equations. The new relation, designated as 
"model PI, " is 


P=x4.]xl.V2lSL.2Q (16) 

L\ I-V L/ 

where 

L = U if |X| s U 
L = |X| if |xl ^ U 

Models Q1 and PI are identical in the principal operating regions. For example, 
in the first quadrant of the P-Q diagram, where P > 0, Q > 0, and X > 0, equations 
(8) and (16) reduce to the following identical relations: 
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Q =X Vl - P 


(17) 


P 



(18) 


Similar results for equations (8) and (16) are obtained by examining the principal regions 
in the second, third, and fourth quadrants of the P-Q diagram. 

Models Q1 and PI have different equations in the center region (|x| ^ U). An 
overlay of Q1 and PI is shown in figure 12. The differences between equations are 
most noticeable for high output pressures (e. g. , | P| > 0. 4). The empirical relation of 
PI may be more realistic than that of Ql, because model Q1 assumes turbulent flow 
through sharp-edge orifices. In practice, leakage passages resulting from spool clear- 
ances do not resemble ideal orifices. The passage surface area, in effect, is large in 
comparison to the flow area. Hence, the more linear ciirves of PI should provide a 
better representation. 

Differences between models are probably insignificant relative to estimated values 
often used for viscous friction, coulomb friction, backlash, and other system param- 
eters. Also, the center region represents only a small portion of the spool's operating 
range. A typical value for the center region would be 3 percent of the maximum spool 
displacement. If duration in the center region is short with respect to the total transient 
time, and if |p| is less than 1, then both models should provide essentially the same 
pressure-flow characteristics. 

Impedance and gain of PI at the origin are derived directly from equation (16). At 
X = 0, P = 0, and Q = 0, we obtain 


3P 2 

Output impedance = — 

3Q U 


3P 2 

Pressure gain = — = — 

ax U 


Flow gain = = 1 

ax 


The small- signal equation at the origin then becomes 


X- Q-ilp = 0 
2 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 
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The same linear small- signal equation results from equation (8) for Ql. There- 
fore, equations (8) and (16) are tangent to each other at the origin. Both models -will 
produce identical system dynamics for small signal disturbances at the origin. 

The appendix gives a computer program for equation (16). Also included in the 
program is a dead zone, or threshold, that simulates the characteristics of overlap 
spools. 


EXAMPLES USING MODEI^ Ql AND PI 

Models Ql and PI were tested with the digital Continuous System Modeling Pro- 
gram (ref. 7) by using either pxire capacitive or pure inertial loads. Pure loads were 
used to simplify the demonstration; they are not valve model requirements. Model Ql 
was tested by using sinusoidal spool motion, a center region (U) of 0. 03, zero overlap, 
and a shunt capacitive load. The input equation was 

X = A sin o)t (23) 

The integral causal relation for the capacitor is 

P=- /*Qdt=-S- (24) 

C •' Cs 

For a linear system, the P and Q waveforms are sinusoidal. Equation (24) then re- 
duces to 


S=ja;C (25) 

P 

Capacitance was set to give the linear case a circular load line on the P-Q diagram. 
This load line is obtained by choosing C - l/a>. Equation (25) then becomes 


^=3 (26) 
P 

According to equation (26), Q and P have equal amplitudes, and Q leads P by a 90° 
phase angle. As a result, the load line on the P-Q diagram is circular and is generated 
by a point that moves in a clockwise direction. 

The load equation obtained from equation (24) for testing Ql then takes the form 
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P = w y* Q dt 


(27) 


Two input amplitudes were used. The smaller signal moved the spool only to the 
edge of the center region. The other signal oscillated the spool to its maximum limits. 
The resultant load lines for steady-state conditions are shown in figure 13 with scales 
that are normalized to amplitude A. 

The curve for the small-amplitude signal (A = 0. 03) is circular with flow and pres- 
sure amplitudes of 0. 03. Deviations from a perfect circle were small. The maximum 
change in the radial direction with respect to a circumscribed circle was only 1. 2 per- 
cent. The circular load line indicates a linear input-output relation. Good linearity for 
small pressure amplitudes is also implied by the characteristics of the curves at the 
Q-axis in figure 12 - they are equally spaced and have equal slopes. 

The curve for A = 1 has an elliptic shape, with its major axis in the second and 
fourth quadrants. This peculiar shape results from variations in flow gain during the 
cycle. Equation (8) indicates that flow gain (9Q/3X) is a function of P. For example, 
if X > U, then flow gain is 1. 413 at P = -1, 1. 000 at P = 0, and 0 at P = +1. Like- 
wise, if the spool displacement is negative and X < -U, then flow gain is 1. 413 at 
P = +1, 1. 000 at P = 0, and 0 at P = -1. Therefore, a large gain change occurs for 
A = 1 when the spool moves through its null position. 

The load curve is generated in a clockwise direction. Consequently, when the 
spool moves through null, for example, the flow gain increases. The result is a sudden 
rate increase in Q. Then, gain gradually decreases as the load curve proceeds through 
the second and first quadrants. The waveforms corresponding to figure 13 are shown in 
figure 14. 

The same two tests were conducted with model PI, except an inertial load was used 
instead of the capacitive load. The load equation was similar to equation (27); that is, 

Q = coy P dt (28) 

The resultant load Unes are shown in figure 15. 

The curve for A = 0. 03 is again circular but evolves in a counter clockwise direc- 
tion. The curve deviates from a circumscribed circle by less than 0. 01 percent, except 
at zero flow, wher-e it changed 0. 6 percent. This transient at Q = 0 is more pronounced 
in the A = 1 cxmve. 

Abrupt changes in pressure occur when load momentum diminishes to zero. Then 
P is unaffected by load and is a function of X only. If X and Q went to zero at the 
same time, P would become zero. Zero pressure, however, did not happen, because 
a small phase lag existed between Q and X. 
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Pressure transients noted in figure 15 did not appear in figure 13. Equation (24) 
explains this behavior. According to equation (24), when Q = 0, then dP/dt is zero. 
Hence, the pressure in an accumulator is constant, as expected, when its flow is zero. 

The waveforms corresponding to figure 15 are displayed in figure 16. Signal dis- 
tortions produced by valve nonlinearities are entirely different from the capacitive load 
curves. A major difference is the shape of the pressvme wave, which reaches an extra- 
high amplitude as the flow approaches zero. High pressure amplitude results from a 
pressure gain (0P/0X) that increases as the load decelerates in the second and fourth 
quadrants of the P-Q diagram. 

A Fourier series of the flow waveform in figure 16 was obtained for the A = 1 
case. The amplitude of the fundamental component is 89 percent of the amplitude ob- 
tained with the linear valve equation (eq. (22)). The phase angle of the fundamental 
component la^ed the linear valve equation by 11°. 

Attempts were made to use model PI with a capacitive load and to use model Q1 
with an inertial load. Solutions for derivative -causality loads could not be obtained 
from the digital Continuous System Modeling Program for either case. This problem 
supports the claim that when the digital computer is used, model Q1 is essential for 
analyzing the effects of fluid compressibility, and model PI is essential for analyzing 
inertial effects. 


CONCLUDING REMARKS 

Programs simulating standard four-way spool valves were developed for two differ- 
ent applications. One was designed for drivii^ capacitive loads, and the other, for 
inertial loads. Both models have demonstrated characteristics that are realistic for 
time-domain simulation of electrohydraulic control systems. 

Model Q1 is a flow- source representation. Input signals are spool position and 
load reaction pressure. Flow is the output signal acting on the load. An underlap 
spool- valve equation defines the leakage through spool clearances. This model is es- 
sentially the same as other published computer programs that have been successfully 
used to simulate valves. Model Q1 should be used in simulations where the capacitive 
effect of entrapped air is the immediate element loading the valve. 

Model PI is a pressure- source representation. Input sigimls are spool position 
and load reaction flow. Pressure is the output signal acting on the load. This model is 
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based on a new constitutive relation for spool leakage. Model PI should be used when 
the immediate load element is either fluid inertia or motor inertia. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 9, 1976, 

778-24. 
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APPENDIX - COMPUTER PROGRAMS FOR MODELING FOUR- WAY SPOOL VALVES 

The following two relations were developed for modeling the pressure-flow charac- 
teristics of four -way spool valves: 

Q = [sign (1 - P)]{0. 5)(X + L) ^\l- p] + [sign (1 + P)](0. 5)(X - L) ^\l + P| (8) 

p =X (2 - mV SM . ^ |i . M\ (16) 

h\ L/l2 hj 

where 

L = U if |x| ^ U 
L = |xl if [x| 2 ; U 

Symbols Q, P, X, and U have normalized values. Variable Q is equal to the 
flow through the output port divided by the maximum output flow at zero differential 
output pressime. Variable P is equal to the differential output pressure divided by the 
maximum differential output pressure that the valve produces for a blocked output line. 
Variable X is normalized spool position, defined as a fraction of maximum displace- 
ment. Symbol U denotes a center operating region of X where spool leakage occurs. 

A t 3 q)ical value for U would be 0. 03; that is, 3 percent of maximum spool displacement. 

Equations (8) and (16) assume that the flow in one output line always equals the flow 
in the other output line; this assumption is generally accepted as a conservative estimate 
for stability analysis. The equations also assume a symmetrical spool, with port areas 
proportional to X. 

Spool stops and overlap characteristics can be included in analog and digital com- 
puter programs. Input signal X is then limited to -1 < X < 1 and modified by a dead 
zone at its center position. The center leakage region remains with overlap spool 
models in order to account for leakage through radial clearances. FORTRAN state- 
ments for equations (8) and (16) that include spool stops and overlap are given in fig- 
ures 17 and 18. Variable V represents spool overlap and is expressed as a fraction of 
maximum spool displacement. Letters D, F, H, K, W, Y and Z are program variables. 
For the zero-overlap case, variable Y is replaced by X. 

Simultaneity of action- reaction signals is an important consideration in the use of 
these models. Bilateral coupling between valve and load requires that the output signal 
and the load reaction signal correspond in time, in order to prevent errors or program 
instabilities. Digital programs produce correlation errors, inasmuch as they can only 
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represent continuous processes by using fixed values during the time-integration inter- 
vals. For accurate results, these intervals must be appreciably smaller than the load- 
time constants. In some cases, this integration interval is so small with respect to the 
overall system response, that the program becomes impractical for use in parametric 
studies. Bilateral coupling presents no problem in analog computers, because simul- 
taneity is a natural product of continuous processes. The valve models are therefore 
more suited for the analog computer. 

The equation preferred for simulating a valve in a control system will depend on the 
load, which may include line dynamics, fluid compliance, actuator leakage, actuator 
inertia, and friction. The valve equation that allows the load equations to have integral 
causality is the one most suited for computer simulations. For example, equation (8) 
is preferred when the effect of entrapped air is simulated as the immediate element 
loading the valve. In this case, P is determined by the fluid accumulated in the cylin- 
ders; that is, the time integral of the difference between Q and the actuator speed. 
Equation (16) is preferred when the immediate load element is fluid inertia in lines con- 
necting the valve to the actuator. Then, P represents the input variable to the load, 
and Q becomes a function of the time integral of P. 

Both valve equations have demonstrated characteristics that are realistic for time- 
domain simulations of electrohydraulic control systems. Each has special applications 
which provide the advantage of using integrators instead of differentiators. 
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Figure 1. - Schematic representation of two-stage electrohydraulic servovalve. 
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Figure 2. - Block diagram of electrohydraulic servovalve. 







Figure 3, - Pressure-flow diagram showing ideal pressure-flow curves of 
critical -center (zero-overlap) four-way spool valves. 
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Figure 10. - Partial block diagram of figure 9 to show program that uses the flow model of the valve. Actuator leakage is zero. 
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Figure 11. - Partial block diagram of figure 9 to show program that uses the pressure model of the valve. Actuator leakage is zero. 
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Y-AMAXl(OwX-V)+AMINl(0.,X+V) 

Z=AMAX1(-1.,AMIN1(1.,Y)) 

F»SIGN(SQRT(ABS(1.-P)),(1.-P))*AMAX1(0.,Z 

H=SIGN(SQRT(ABS(i.+P)),(l.+P))*AMINl(0.,Z 

Q=F+H 


,0.5*(Z+U)) 

,0.5*(Z-U)) 


Figure 17. - Block diagram and FORTRAN statements of equation (8). 
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